Inhomogeneous High Harmonic Generation in Krypton Clusters 
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High order harmonic generation from clusters is a controversial topic: conflicting theories exist, 
with different explanations for similar experimental observations. From an experimental point of 
view, separating the contributions from monomers and clusters is challenging. By performing a 
spectrally and spatially resolved study in a controlled mixture of clusters and monomers, we are able 
to isolate a region of the spectrum where the emission purely originates from clusters. Surpringly, the 
emission from clusters is depolarized, which is the signature of statistical inhomogeneous emission 
from a low-density sources. The harmonic response to laser ellipticity shows that this generation is 
produced by a new recollisional mechanism, which opens the way to future theoretical studies. 

PACS numbers: 



High harmonic generation (HHG) f]\ refers to the in- 
teraction of high intensity laser light with matter, which 
leads to the emission of broadband coherent radiation 
in the extreme ultraviolet domain. HHG from clusters 
is considered a promising light source, showing higher 
emission frequencies [J. HHG is also a spectroscopic 
tool to extract structural and dynamical information on 
the emitting medium from the properties of the harmonic 
radiation (spectrum, phase and polarization state). This 
technique, which has been used to probe atoms [^,'5'] and 
small molecules [6] , relies on the basic mechanism of HHG 
described as a three step model 0,0,0: First, a bound 
electron escapes in the strong laser field through tunnel 
ionization. Second, the electron is driven away then ac- 
celerated back towards the parent ion. Finally the elec- 
tron recombines radiatively with the parent ion. This 
recombination encodes the structure of the medium in 
the emitted light. The extension of this technique to the 
case of clusters would allow investigation of strong field 
processes in many-body systems, the role of multielec- 
tron effects and monitoring of cluster dynamics through 
the harmonic signal 0, [lo| . 

The exact mechanism of HHG from clusters is still de- 
bated. Various extensions of the three step model have 
been proposed |lll-ll4|. The dominant channel gener- 



ally considered is ionization and recombination to the 
same atom (atom-to- itself). Since clusters are dense me- 
dia, there is also a possibility of recombination to neigh- 
bouring; ions [13, [la. This atom-to -neighbour emission 
can produce incoherent radiation due to a lack of phase- 
locking between the two atomic wavefunctions [14]. An- 
other contribution to harmonic emission may come from 
a wavefunction partially delocalized over the whole clus- 
ter, from which electrons tunnel out of and to which 
they recombine coherently (cluster-to-itself). In addition 
to these recollisional mechanisms, very different mech- 
anisms could also co-exist as is the case in overdense 
plasmas [li^ or bulk crystals [16]. From an experimen- 



tal point of view a particular difficulty consists in dis- 
entangling the harmonics produced by different species 
(monomers and different size clusters) and possibly dif- 
ferent mechanisms. 

In this letter we present a detailed experimental study 
of HHG from clusters which aims at disentangling the 
harmonic signal from clusters and identifying the mecha- 
nisms at play. By studying the spectral and spatial profile 
of the harmonics as a function of pre-expansion gas tem- 
perature, we identify one region corresponding to pure 
emission from clusters. A simulation verifying the re- 
sults from a polarization measurement confirm that only 
few emitters contribute to this region, which is consistent 
with the low density of large clusters in the generating 
medium. The fast decay of the harmonic signal with laser 
ellipticity indicates that these harmonics are produced by 
a recollision mechanism which strongly depends on the 
cluster size, as expected in a cluster-to-itself picture. 

We generate high harmonics in a supersonic gas jet, 
employing a valve (Even-Lavie [17]) pulsed at IkHz with 
a d = 150 /im diameter trumpet nozzle and a jet expan- 
sion half angle of a ^ 10°. The backing pressure was 
kept constant at po = 20 bar, while the pre-expansion 
nozzle temperature Tq varied. Cluster formation can 
be described by the empirical Hagena scaling parameter 

T* = k ^^^^^2^1 Po [18], with a condensation parame- 



ter k equal to 2890 for krypton [18]. The average cluster 
size can be estimated by [l9|: N = 33 (rVlOOO)^-^^ In 
our experiment the cluster size was varied by changing 
the nozzle temperature between 350 K {N ^ 33000) and 
510 K (TV ^ 4000). The standard deviation SN, of the 
cluster size distribution is broad and leads to significant 
contribution of small cluster, which can be estimated by 
SN/N = 45% (l9|. 

We used the Aurore Ti:Sa-Laser system from CELIA 
which delivers 7 mJ, 35 fs, 800 nm pulses at 1 kHz. The 
laser was focused using a 37.5 cm focal spherical mirror 
at 4.8 mm from the nozzle exit, which ensures proper jet 
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shell ionzation [24*] and leaving the freed electrons in a 
quasy bound cloud behind, can be responsible for the 
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FIG. 1: (a) High harmonic spectrum at 350 K. The vertical 
lines correspond to electronic transitions in krypton ions, (b) 
Normalized difference signal of a spectrum at 350 K and a 
spectrum at 510 K. The signal from the pink areas dominates 
at 350 K. 



thermalization. The harmonic spectrum was dispersed 
using an abberation-corrected concave gold grating. The 
imaging detector consisted of dual microchannel plates 
coupled to a phosphor screen and a 12 bit CCD-camera. 
Each image was averaged over 15000 laser shots. 

Figure [T]^a) shows a harmonic spectrum obtained at 
a laser intensity of IxlO^^W/cm^. The horizontal line 
around zero divergence angle is due to scattered light 
from lower harmonics. In addition to the harmonic spec- 
trum, we observe isotropic radiation as vertical lines, cor- 
responding to plasma lines [2o|, [Hj , which originate from 
electronic transitions in excited krypton ions. Plasma 
lines at 38.7 nm and lower wavelengths correspond to 
electronic transitions in at least five times ionized kryp- 
ton ions [22]. A strong line observed at 43.4 nm refers 
to an electronic transition in Kr^+ [23|. The laser inten- 
sity employed is not high enough to field ionize krypton 
atoms five (1^=78.5 eV) or seven times (1^=125.8 eV). In 
turn, impact ionization within clusters leading to inner 
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formation of such highly ionised Kr 

Apart from plasma lines, the spectrum from Fig. [T]^a) 
does not show specific features that could be the signa- 
ture of the presence of clusters in the medium. In order to 
identify the contribution of clusters we performed a dif- 
ferential measurement: we monitored the spatially and 
spectrally resolved evolution of the signal as the nozzle 
temperature increased from 350 K to 510 K, which de- 
creases the cluster average size. Fig. [IJb) reveals clear 
regions with distinct behaviors: a low energy compo- 
nent (up to harmonic 25) which increases with temper- 
ature, and an off- axis component (dominant from har- 
monic 23 to 31) which presents a higher cutoff and de- 
creases monotonously with temperature (Fig.fljb)) as the 
plasma line. In the following we will focus on the latter 
component, whose behavior is correlated to that of the 
plasma lines and thus to the presence of large clusters. 

In order to characterize the different mechanisms that 
may be at play in the harmonic emission, we performed a 
polarimetry measurement [26|. We used an unprotected 
silver mirror in combination with a gold grating as a fixed 
polarizer and rotated the laser polarization with a zero- 
order half waveplate. The signal of each pixel of our 
detector is monitored as a function of the angle of the 
half waveplate resulting in a Mains' law with a cos'^ de- 
pendence. After subtracting the background, a Fourier 
transform is performed to extract the amplitude of the 
oscillatory component, which is then normalized to the 
sum of the Fourier transform. This procedure provides 
the amplitude of the oscillation and therewith the de- 
gree of linear polarization. Fig. [2fa) depicts the spatio- 
spectrally resolved degree of linear polarization obtained 
at 350 K. As expected, the plasma lines are clearly unpo- 
larized. More surprisingly, the same region we attributed 
to cluster emission in Fig. [Hb) shows a remarkably low 
degree of linear polarization. A low degree of linear polar- 
ization is very unexpected for high harmonic generation 
by a linear laser field [26]. Note that we systematically 
observed this low degree of linear polarization in different 
conditions and various species (CO2, argon and mixtures 
of CO2 /krypton), but not in the absence of clusters using 
an effusive gas jet. 

The polarization direction of the high harmonics is set 
by the recollision direction of the electrons and the elec- 
tronic structure of the ground state to which these elec- 
trons recombine. Within the strong field approximation, 
in linear polarization the recollision direction is parallel 
to the laser field. In a centrosymmetric medium the har- 
monics are thus necessarily polarized along this direction. 
To go beyond the strong field approximation, we have 
performed Classical Trajectory Monte Carlo (CTMC) [H 
calculation of the recollision angles in atoms at intensi- 
ties of 1.2 xlO^^W/cm^. These calculations take into 
account the influence of the ionic core on the electron 
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FIG. 2: (a) Polarimetry map depicting the oscillation am- 
plitude of the the Mains' law for each pixel at 350 K. (b) 
(b) Simulated harmonic spatial profile for a single laser shot 
(left), 1000 shots (middle), and degree of linear polarization 
for 1000 shots (right). The first row corresponds to HHG in 
atoms and the second row to HHG in large clusters. 

trajectories, which broadens the distribution of recolli- 
sion angles. After summing the recolliding electrons over 
the polar angle, we obtain a recollision angle distribu- 
tion peaked around ±6° [HI, [l^. In clusters, given the 
high charge states observed in our experiments, we ex- 
pect the ionic potential to play a more important role 
[28|. Futhermore screening effects can inhomogenize the 
electric field over the cluster [29]. These two effects will 
lead to a broader recollision angle distribution. 

For understanding the origin of the apparent depo- 
larization of harmonic emission, we performed simple 
simulations of the macroscopic generation process as- 
suming an infinitely thin medium. We randomly dis- 
tribute N emitters on a square grid of 100*100/im which 
represents the generating medium. Each emitter lo- 
cated in (x^y) radiates an electric field Eq{x^y) = 
/o(x,7/)^^///2e^^'2^o(^'^)^r(x,?/), where Io{x,y) is the fun- 
damental intensity distribution, QEff the effective non- 
linearity of the harmonic emission (typically 5), aq the 
intensity dependence of the harmonic phase (a ~ 10 x 
10"-^^ cm^/W for the end of the plateau), and u{x^y) 
is unit vector with a direction randomly picked in the 
distribution of recolliding angles. For clusters the lat- 
ter is assumed to be three times broader than the one 



obtained for atoms by using CTMC. We calculate the 
electric field resulting from the coherent sum of the con- 
tributions from the N emitters in the far field by Fourier 
transforming the near field profile. Two extreme cases 
are considered: emission from an ensemble of clusters 
and from an ensemble of monomers with a density of 
10^^ cm~^. Assuming that only one atom out of 10^ 
emits high harmonics (due to the recombination prob- 
ability) and a medium thickness of 100 /im, we get 10 
atoms to distribute on our grid. For clusters, the den- 
sity is much lower: assuming a medium with 80% of the 
atoms forming clusters of 30000 monomers, the density 
is ~ 3 X 10^^ cm~^ [19], which leads to a few hundred to 
thousand emitters with an emission probability the same 
or higher than that for atoms. In the following we will 
consider 1000 emitters and we checked that the results 
were qualitatively similar using 100 to 10000 emitters. 

These simulations show that while in the case of atomic 
emission the harmonic far field profile is well defined in a 
single laser shot, the profile obtained from clusters is very 
inhomogeneous and fluctuates significantly from shot to 
shot because there are too few emitters to obtain a nice 
constructive interference on axis and destructive off axis 
(Fig. Efb)). The polarization state from clusters (not 
shown) is also very inhomogeneous, showing important 
polarization angles and ellipticities. When averaging over 
1000 shots, the spatial intensity profile becomes much 
smoother and the polarization angle is quasi homoge- 
neously equal to zero. We also calculate the degree of 
linear polarization of the resulting light, by defining a 
Malus law out of the Stokes parameters and extracting its 
normalized oscillatory component. While the polariza- 
tion appears perfectly linear in monomers, the emission 
from clusters appears depolarized. These results confirm 
that the unpolarized cutoff harmonics (Fig. [IJb)) are 
a fingerprint of HHG from large clusters, which have a 
low density in the jet and thus produce inhomogeneous 
harmonic emission. 

In order to further elucidate clusters as a class of emit- 
ters, we studied the sensitivity of the high harmonic 
signal to the ellipticity e of the generating laser field. 
The main polarization axis was kept vertical by employ- 
ing a fixed zero-order quarter waveplate behind an ad- 
justable zero-order half waveplate. The ellipticity was 
varied from -0.45 to 0.45 recording 49 points. As e in- 
creases, we observe an exponential decay of the whole 
harmonic spectrum, which is consistent with a recolli- 
sion picture of HHG. The decay rate /3 of the harmonic 
intensity / with ellipticity can be evaluated by a Gaussian 
fit: I (X exp{—j3e^). We extract this decay rate for each 
pixel of the detector and obtain the results shown in Fig. 
[31(a). As in Fig. [Hb)) the same areas can be identified, 
the off-axis cutoff with /3 > 45 and on-axis plateau with 
/3 < 45. Remarkably the nozzle temperature increases 
the decay rate in the first case, while it hardly affects 
j3 in the second case (Fig. [31(b)). When extracting the 
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FIG. 3: (a) Harmonic decay rate p at 510 K. (b) Evolution of 
P with temperature. The values are an average over the areas 
depicted above. 



slope for each pixel, a positive slope is observed only for 
the off-axis cutoff region. 

The exact mechanism of high harmonic generation in 
clusters, responsible for the emission of the highest har- 
monics, has yet to be determined. First, let us mention 
that the laser intensity and electron density are too low 
to produce significant coherent wake emission from the 
plama [iSj. Second, the fast decay of the harmonic signal 
with ellipticity indicates that the process is recollisional 
and is thus different from what was recently reported in 
the case of bulk crystals [16], which still shows a signif- 
icant signal at e = 0.5. Third, recombination to neigh- 
bouring ions [l^ can also be excluded as it would lead 
to lower /3- values. We therefore suggest another mecha- 
nism of HHG in clusters, namely tunnel ionization from 
a partly delocalized electron wavefunction and recombi- 
nation to this wavefunction (cluster-to- itself). Even lo- 
calized electronic states in a Van der Waals clusters [30| , 
could be driven by the strong laser field from one site 
to another, ending up in a partially delocalized wave- 
function after a few optical cycles jsij. This process 
is accompanied by a significant amount of ionization of 
the cluster (as observed from the plasma lines) which 
will increase its ionization potential and consequently the 
higher cutoff observed for the cluster emission. Record- 
ing the harmonic signal versus e, one obtains the cross- 
correlation between the recolliding wavefunction and the 
ground wavefunction. In the "cluster-to- itself picture 
the electrons of interest are expected to tunnel mostly 
from the surface of the cluster, so that the initial width of 
the electron wavepacket after tunneling is proportional to 
the cluster size. The wavepacket spreads laterally during 



acceleration outside the cluster - the smaller the initial 
wavepacket, the stronger the spread. The decay rate, /3, 
is determined by the width of the recolliding wavepacket 
and the extension of the initial wavefunction, i.e. the 
cluster size. For large enough clusters the latter will be 
dominant, so that the decay rate is expected to decrease 
with increasing cluster size. This is consistent with the 
observed temperature dependence of H27 in Fig. [3fb). 

In conclusion, by performing a 2D spectro-spatial anal- 
ysis, we are able to disentangle several contributions 
to the harmonic signal from a mixture of clusters and 
monomers. The high-energy off- axis emission shows 
three clear features: a higher cutoff of the harmonic emis- 
sion, a higher decay rate with ellipticity when the average 
cluster size is reduced, and a very low degree of linear po- 
larization. We attribute these effects to a high harmonic 
generation mechanism in which delocalized electrons tun- 
nel from a cluster and recollide coherently to the whole 
cluster. This hypothesis needs to be further investigated 
using the appropriate theoretical tools. This work could 
be extended by performing a complete polarimetry study 
of the harmonic emission in elliptical laser fields, which 
was recently shown to be a sensitive probe of the in- 
fluence of the ionic potential in HHG [27|. This would 
enable us to determine the process responsible for the 
generation of the lowest harmonics and possibly further 
differentiate the contribution from monomers and from 
atom- to- atom emission in clusters. The use of advanced 
2D experimental characterization to disentangle contri- 
butions from different classes of emitters and different 
mechanisms can be applied to many situations in high 
harmonic spectroscopy, and should in particular be use- 
ful for other cases of HHG in inhomogeneous media, like 
metal plasmas [32]. 
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